Driving linear systems to non-classical states with the help of noise 
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We study the possibility of taking bosonic systems subject to quadratic Hamiltonians and a noisy thermal 
environment to non-classical stationary states by weak Gaussian measurements and conditioned linear driving. 
We derive general analytical upper bounds for the single mode squeezing and multimode entanglement at steady 
state, depending only on the Hamiltonian parameters and on the number of thermal excitations of the bath. Our 
findings show that, rather surprisingly, larger number of thermal excitations in the bath allow for larger steady- 
state squeezing and entanglement if the efficiency of the Gaussian measurements conditioning the feedback loop 
is high enough. Such efficiencies are included in our exact treatment, which allows us to determine efficiency 
thresholds for the noise-enhancement of quantum resources to take place. 
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All quantum technologies hinge on establishing controlled 
interactions between different constituents of quantum sys- 
tems whilst reducing unwanted interactions with an environ- 
ment, which give rise to decoherence. In dealing with environ- 
mental decoherence, two main paradigms have emerged over 
the last fifteen years: one may either attempt to decouple the 
relevant, logical degrees of freedom from the environment by 
various techniques (e.g., decoherence free subspaces yj], error 
correction yj], dynamical decoupling yD), and then proceed 
to process the quantum information coherently (e.g., in gate- 
based models of quantum computation, through unitary oper- 
ations), or one may try to manipulate the noisy, non-unitary 
evolution of the system directly, tailoring it to suit one's aims. 

The second viewpoint, which one might broadly refer to 
as the 'dissipative' approach to quantum information process- 
ing, has a long tradition, going back to the first proposals for 
reservoir engineering [|4|, and has recently been compounded 
by the design of a model for dissipative, non-unitary quan- 
tum computation ||5|]. It has hence been repeatedly shown, 
in various contexts and settings, that working with the en- 
vironment rather than against it may lead to forms of coop- 
eration whereby the environment contributes to enhance cer- 
tain coherent tasks performed on the system, often in a rather 
counterintuitive manner ll6l425ll . This Letter is the account of 
such cooperation in the broad setting of controlled dissipative 
dynamics in linear Gaussian systems 111 127143311 . We shall 
consider a system of n bosonic modes subject to a quadratic 
Hamiltonian and to dissipation in a thermal environment with 
average excitation number N, and show that the maximal 
squeezing and entanglement achievable by continuous linear 
feedback control grows with N, that is with the temperature 
of the bath. We will apply our results to various quadratic 
Hamiltonians and also study quantitatively the role played by 
the efficiency of the weak measurements that condition the 
feedback loop. 

Notation and background. - We consider a system of n 
bosonic modes described by the vector of canonical opera- 
tors R = {xi,pi, . . . , i;„,p„)^, with commutation relations 



encoded by the anti-symmetric symplectic form il, as per 
[Rj, i?fc] = iOjfe (?i = 1 throughout the paper). 

Being comprised of Gaussian noise, Gaussian averages, 
and Hamiltonian evolutions of the first ('linear driving') and 
second ('canonical', or 'symplectic') order in the canonical 
operators, our dynamics will only involve Gaussian states, 
which are entirely described by first and second statistical 
moments of the canonical operators ll34ll . The second mo- 
ments of a Gaussian state g, in particular, will be repre- 
sented by a 2n X 2n 'covariance matrix' (CM) cr: ajk ~ 
Tii{Rj,Rk}g) - 2Ti{Rjg)Tr{Rkg), which satisfies the well 
known Robertson-Schrodinger uncertainty relation: 



iO > 



(1) 



This is a necessary and sufficient condition for a CM to repre- 
sent a physical Gaussian state g 135 1. 

Let us begin by considering the most general time- 
independent quadratic Hamiltonian acting on the system: 

H = (l/2)RTi7R , 

where the 'Hamiltonian matrix' H isa generic symmetric ma- 
trix. We will later on modify the Hamiltonian to include a 
time-dependent linear term which will exert the feedback ac- 
tion on the system. Besides, the system is linearly coupled 
to a large uncorrected Markovian bath, whose action can be 
described, in the standard input-output formaUsm dj] |36|], by 
a stochastic Langevin equation for the canonical operators R: 
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where A is the 'drift matrix', given by A = (riH — k12„)/2 
(111 standing for the d-dimensional identity matiix), k is the 
loss rate of the system and Ri„ is the time-dependent vector 
of the bath's operators coupled to all the different system's 
modes [IJ. We will consider the most general case of Gaus- 
sian white noise, with bath's correlation functions encoded in 
the diffusion matrix D according to K{Rin,j{t)Rin,k{t')) = 
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FIG. 1. Heuristic pliase-space representation of an optimal linear 
feedback action. The unconditional state is a Gaussian average, with 
CM <T, of conditional Gaussian states with the same CM cTc and 
different centres in phase space (a). The optimal Markovian choice 
for the linear driving term, represented by gray arrows in (a) and by 
u{t) in the Hamiltonian, cancels the first moments of the conditional 
state, thus making it coincide with the unconditional averaged one 
(b). 



Djk5{t—t'). Finite temperatures of the environment are incor- 
porated in our treatment by choosing D = ® i=i(l + 2A'j)l2, 
where Nj represents the number of the thermal excitations in 
the bath of mode j OTIl . The free (in the sense that no moni- 
toring or feedback actions have been introduced yet) dynamics 
of the second moments under such conditions is described by 
da/dt = Act + crA^ + D. If the system is stable, in the sense 
of admitting a steady state, it must be [A + A^) < 0, which 
we will assume in what follows. Also, from now on we will 
set K = 1, and use the loss rate as a unit of frequency. No- 
tice that, to simplify our notation, we have assumed the same 
loss rate k affecting all modes. This condition could be re- 
laxed and, in fact, all our general results will be presented in 
terms of generic drift and diffusion matrices, which allow for 
generic loss rates. 



Let us now define the output field Rout, by the boundary 
condition Rout = ^/kR — Rm llll. Monitoring the output 
field Rout, which is effectively coupled to the system's field, 
corresponds to performing a weak measurement on the latter: 
we will now include in our dynamics the most general Gaus- 
sian "general-dyne detection" on the field Rout llll, encom- 
passing all homodyne detections, both direct and resorting to 
ancillary modes (and hence heterodyne detection too). Such 
measurements are the most general allowing for a continu- 
ous, though stochastically fluctuating, monitored evolution of 
the system if one assumes, as customary in the framework of 
'continuous' feedback control, that the monitoring of the envi- 
ronment happens on time-scales which are much shorter than 
the typical system's response time. The detailed description of 
an arbitrary general-dyne measurement and of the associated 
monitored dynamics is rather cumbersome. Fortunately, how- 
ever, our results can be presented by introducing only minimal 
technicalities (see the supplemental material for a systematic 
treatment). It will suffice here to say that the effect of any 
general-dyne measurement of Rout can be described in terms 
of an unravelling matrix U. The knowledge of U allows one 



to reconstruct the specific Gaussian measurement by which 
the environment is being monitored. Furthermore, we allow 
for a dependence of U on the parameter rj, quantifying the 
overall efficiency of the measurements (in the sense that, be- 
fore any measurement, a proportion (1 — 77) of the signal, with 
< 77 < 1, is lost). We will hence be able to include efficien- 
cies in our analytic treatment. 

The continuous monitoring of the output field is recorded 
in the current y(t), and leads to a conditional dynamics 
which does not alter the Gaussian character of the system's 
state (here, "conditional" refers to the conditioning due to the 
knowledge of the weak measurements' outcomes). In the fol- 
lowing, we shall distinguish between the conditional state of 
the system Qc, with CM cTc, and the time-averaged, 'uncondi- 
tional' state g = (1/At) J^ gc{s)ds, with CM er, where 
At is an integration interval much larger than the typical time- 
scale of the stochastic fluctuations of the measured current 
13811 . As we mentioned, the dynamics of the conditional state 
Qc is Gaussian, with stochastic fluctuations (depending on the 
measured current) affecting the first moments, but an entirely 
deterministic evolution for the matrix of second moments cTc- 
This fact, as we will see, is essential to our discussion. In fact, 
the white-noise fluctuations of the first moments are so fast 
that one is left with the average, unconditional evolution of 
the quantum state to all practical purposes. But, as depicted in 
Fig.[T] the unconditional state is just a Gaussian state resulting 
from the average of conditional Gaussian states with the same 
CM <Tc and different first moments (centres of their positions 
in phase space). It is very easy to see that, under such an av- 
erage, all the figures of merit we are going to consider can 
only decrease. Hence, for given general-dyne measurement, 
the best case scenario for any of our figures of merit would 
be one where the fluctuations of the first moments cancel out 
and the average unconditional state coincides with the condi- 
tional state. It turns out that such a situation can always be 
arranged by adding a linear Markovian feedback action to the 
Hamiltonian: 



Hf = -BJnBY{t) 



(2) 



where y(t) is the general-dyne current and i? is a matrix com- 
pletely determined by the unravelling matrix U (see the sup- 
plemental material for the explicit expression of B). Marko- 
vian feedback is therefore always optimal to our aims and we 
will hence restrict to it in the following. Before proceeding, let 
us briefly mention that the dynamics of the averaged, uncon- 
ditional second moments under a linear Markovian feedback 
action like that of Eq. ( fTSI l can still be treated analytically and 
is of the form da/dt = A'ct + aA''^ + D' (the modified drift 
and diffusion matrices are given in the supplemental material). 
In view of the above, in order to optimise the steady state 
squeezing or entanglement, one has just to optimise the rele- 
vant figure of merit for the conditional state Qc, and then ap- 
ply the Markovian feedback strategy that ensures g = Qc (see 
Fig. [Til. The optimization over the set of conditional states 
does not need to go into the details of the conditional dynam- 



ics but can instead be tackled by resorting to a general math- 
ematical result: given drift matrix A and diffusion matrix D, 
a CM <Tc is a stabilising solution of the deterministic condi- 
tional dynamics of the second moments if and only if 12911 



Ac 



CTrA^ + D > 



(3) 



General results. - We first present two very general bounds 
following from Inequalities (dJ and ^ and from the optimal- 
ity of Markovian feedback discussed above and illustrated in 
Fig. [T] (see the supplemental material for the proofs). In the 
following, oi'- will stand for the j-th smallest eigenvalue of 

the positve definite matrix — (^ + A^), while 5^ will stand for 
the j-th largest eigenvalue of the positive definite matrix D. 
We start by considering the achievable steady-state squeezing, 
defined as the smallest eigenvalue A{ of the CM cr. 

Proposition 1 (maximal squeezing). - Let a be the CM 
of a steady-state achievable by continuous weak general-dyne 
measurements and linear driving in a system of bosonic modes 
subject to a drift matrix A and Gaussian white noise with a 
diffusion matrix D. The squeezing A| is bounded by 



5] 



(4) 



Next, we turn our attention to the logarithmic negativity 
Ej^, an entanglement monotone which provides one with an 
upper bound to the distillable entanglement, can be evaluated 
analytically for Gaussian states, and is a good quantifier of 
quantum correlations for systems which allow for no bound 
entanglement (like Gaussian bipartitions of 1 versus (n — 1) 
modes or bisvmmetric |39ll Gaussian bipartitions of any num- 
ber of modes) ^ EH W^ 

Proposition 2 (maximal entanglement). - Let g be the CM 
of a steady-state achievable by continuous weak general-dyne 
measurements and linear driving in a system of bosonic modes 
subject to a drift matrix A and Gaussian white noise with 
a diffusion matrix D. The logarithmic negativity Ej\f{Q) of 
any 1 versus (n — 1) modes or bisymmetric bipartition of g is 
bounded by 



Ej\f{g) < max 



0,log2 



(5l + 4 



,2\/ala|, 



(5) 



A remarkable feature of both our bounds is that they in- 
crease (somewhat loosely, we will refer to Aj getting smaller 
as an 'increase' in the squeezing) if the largest eigenvalues of 
the diffusion matrix D increase, which characterises a noisier 
environment. As mentioned above, if one considers a simple 
thermal environment, one has 5l — 82 — '^ + 2N^, where N^ 
is the largest number of thermal excitations in an environmen- 
tal degree of freedom. As we shall see, our bounds are actu- 
ally tight in several important cases, where they represent the 
actual maximal values achievable. Hence, our findings show 
that the maximal achievable entanglement increases with the 
temperature of the bath. 



p 




P 




f 


X 


— ^j\ h -^ X 


I 


J 


■^ ^ 




(a) ^ 


-^!h. 


(b) ■ •• 


n 



FIG. 2. Heuristic phase space picture of the noise-enhancement of 
the optimal feedback action. The feedback squashes the thermally 
broadened unconditional steady state CM CTth, turning it into the 
squeezed CM a. By Heisenberg principle, the squashing is limited 
by the inverse of the thermal uncertainty in the orthogonal quadra- 
ture, which increases with increasing noise. 



This apparently counterintuitive property can be illustrated 
and understood by considering the feedback action on the 
squeezing of the unconditional state of a free single bosonic 
mode 14311 . As we will see, in this case the optimal procedure 
to obtain squeezing consists in monitoring the environment 
along a given phase space direction and then in systematically 
driving the expectation value of the monitored quadrature to 
zero. As illustrated in Fig. |2] this produces an uncertainty 
contraction for that quadrature, while the conjugate, orthog- 
onal quadrature is entirely unaffected. Hence, by Heisenberg 
principle, the achievable squeezing is ultimately limited by the 
inverse of the uncertainty in the orthogonal quadrature, which 
clearly increases with the available thermal energy of the bath. 
In a sense, this is a case of reservoir engineering where the ef- 
fect of the bath is 'squashed' 14411 . rather than squeezed, by 
means of continuous measurements. 

If one is interested in optimal squeezing, this thermal en- 
hancement can be obtained by measuring and acting locally 
on a single quadrature, while the generation of optimal en- 
tanglement will generally require nonlocal measurements. It 
should however be noted here that linear feedback does allow 
for an increase in steady-state unconditioned entanglement 
even with local measurements, if the Hamiltonian couplings 
between the modes are strong enough 113 311 . Let us also men- 
tion that we have determined precise necessary and sufficient 
conditions for the bounds to be achievable (these are included 
in the supplemental material), and that, whenever the bound is 
achievable, the optimal steady-state is pure, because the satu- 
ration of the uncertainty relation ([!} is implied. In such cases, 
the optimal feedback strategy not only maximises a figure of 
merit but also stabilises a pure state, regardless of how noisy 
the environment may be. 

In what follows, we shall focus on specific examples and 
discuss findings of direct applicative interest. The reader may 
refer to the supplemental material for further technical details, 
such as explicit expressions for achievable optimal states and 
optimal unravellings. 

Free systems. - Let us start with the simple case where 
no Hamiltonian is present {H = 0, which in practice corre- 
sponds to considering a system in the rotating frame, and to 



having all the measurements' phase references rotate accord- 
ingly). Henceforth, we will always set Si = S2 = i- + 2A^{' 
(phase-insensitive thermal noise). Without any feedback ac- 
tion, the steady-state clearly corresponds to a thermal state 
without squeezing nor entanglement. On the other hand, the 
bounds on the squeezing and logarithmic negativity achiev- 
able via feedback read, respectively, A| > 1/(1 + 2A^j) and 

E^ < log2 1 1 + 2Ni J . The bound on the squeezing is al- 
ways achievable by a simple feedback strategy; one has to 
perform a weak measurement on a certain (arbitrary) quadra- 
ture Xij, interacting with the noisiest bath, and then drive the 
orthogonal quadrature with the corresponding current. The 
resulting state will be optimally squeezed in the measured 
quadrature i^. A very instructive and consequential curios- 
ity arises here: the efficiency of the optimal unravelling is not 
77 = l,butrather7? = N^{l+N^)/[1+Nf{l+N^)],v^hichis 
always smaller than 1 for N^ > 0. This is a very explicit case 
of information-disturbance trade-off, whereby a weaker mea- 
surement (as is the case for lower efficiency) holds the virtue 
of causing less disturbance on the state, and might be desirable 
at times. This fact, besides its fundamental interest, might also 
have practical consequences, in systems where very few ther- 
mal excitations are the dominant source of noise, as could be 
the case at terahertz frequencies in solid-state and optical sys- 
tems at room temperature 11451 14611 : if N^ ~ 1 (corresponding 
to about 30 THz), then the optimal detection efficiency is just 
r] ~ 0.66. It should however be noted that the optimal effi- 
ciency raises very quickly to 1, being already ?/ ~ 0.86 for 
Nf ~ 2 (around 15 THz). 

As for the entanglement, we can show that in the two- 
mode (n = 2) case, if the two baths have the same tem- 
perature N^, the bound can be saturated. One of the possi- 
ble optimal unravellings and Markovian strategies correspond 
to driving the quadratures {xi + X2) and {pi — P2) with 
the currents obtained by monitoring the conjugated quadra- 
tures {pi + P2) and {xi — X2), respectively. In this case 
an ideal homodyne measurement with unit efficiency is re- 
quired to achieve the bound, and we hence considered explic- 
itly also the case where this strategy is applied with inefficient 
measurements. Then, the logarithmic negativity achieved 

for efficiency 77 is given by £V = log2 ( 1 + 2A^J' j — 
log2 ( 1 + 4iVj'-(l - 7?) + 4:Nf{l - 77) ] . By inspecting this 
equation one observes that, for a given temperature iV|, one 



the case of degenerate parametric down conversion which can 
be described, in interaction picture, by the quadratic Hamilto- 
nian H = \ix1P2 + P1X2) between two modes at the same 



2(T+]vjy 



such that entan- 



can define a threshold value 77th 
glement is obtained only for efficiencies 77 > 77th. We no- 
tice that rjth is always greater than 1/2 and monotonically in- 
creases with temperature towards the maximum value corre- 
sponding to a perfect homodyne measurement. 

In the most general case, when N{ 7^ iVj , the bound cannot 
always be saturated. However, we were able to find an unrav- 
elling such that the steady state is a pure two-mode squeezed 

state with logarithmic negativity £V = log2 ( 1 + 27V2 
Parametric Hamiltonians. - We now move on to consider 



frequency 114711 . such that the average number of thermal ex- 
citations in the two modes are the same and set equal to N: 
Z? = (1 + 2iV)l4. We shall impose stability by bound- 
ing the interaction strength: x < 1/2. This set of dynami- 
cal parameters allows for the perfect saturation of the bound 
on the entanglement, and thus for the analytical optimisation 
of the achievable logarithmic negativity, which is given by 
Ej\f ~ log2(l + 2A^) — log2(l — 2x), to be compared with the 
free steady state value E^j^^ = log2(l + 2x) - log2(l + 2N) 
(what would be obtained in the absence of monitoring and 
feedback action). This is possibly the most apparent exam- 
ple of noise-enhanced performance in our study: while the 
free steady state logarithmic negativity decreases with the 
temperature, as one should expect, its optimised counterpart 
increases with N. The closed-loop control we considered 
is capable of retrieving information from the output channel 
and turning the phase insensitive thermal energy into correla- 
tions between the modes. The optimal feedback strategy has 
been exactly determined in this case as well, and corresponds, 
again, to driving the quadratures {xi + X2 ) and {pi — p2 ) with 
the currents obtained by monitoring the conjugated squeezed 
quadratures {pi + P2) and {xi — X2), respectively. In this 
case too a perfect homodyne measurement is required and 
one should hence consider the effect of the efficiency 77 on 
the achievable optimal entanglement. The conditions on the 
measurement efficiency for the feedback loop to be able to im- 
prove the generation of entanglement are rather strict, and be- 
come steeper and steeper as the noise increases. For x = 0.3 
and A^ = 1, where the steady state in absence of feedback 
is unentangled, 77 > 0.8 is needed to generate any entangle- 
ment between the two set of modes. This threshold increases 
to 0.92 for N = 2.5 and to 0.98 for N = 10. As already men- 
tioned above, these are hence the typical values of excitations 
where linear feedback control might really make a difference 
in the generation of pure entangled states of continuous vari- 
able systems. 

Position couplings. -Let us also consider the other paradig- 
matic quadratic Hamiltonian H ~ x 2;ii;2 where, again, the 
two modes have the same frequency and the average num- 
ber of thermal excitations in the two modes is set equal to 
N. This is the typical form of linearised interaction involv- 
ing material harmonic oscillators, such as Coulomb-coupled 
trapped ions or mirrors coupled to cavity field modes in opto- 
mechanical set-ups. In this case, one observes that the bound 
obtained by means of our theorem is not achievable, since the 
conditions outlined in the supplemental material are not sat- 
isfied. However, we were able to identify a family of two- 
mode squeezed states which are stabilising solution of the 
conditional dynamics and whose logarithmic negativity reads 
Ef/ = log2(2(l + 2N)) — log2(x + 2). As is apparent from 
the formula, the entanglement of these states increases with 
the thermal excitations number N, which proves that a ther- 
mal environment can boost the entanglement of xx coupled 



oscillators as well. 

Conclusions. - In this Letter, we have shown that linear 
feedback may be employed not only to stave off the effect 
of thermal noise and achieve pure steady states ll48ll49ll . but 
also to convert the thermal energy of a bath into enhanced 
squeezed or entangled resources, of paramount importance to 
precision measurements ll50t - l52ll . quantum information pro- 
cessing 15 3115411 . and quantum communication with continu- 
ous variable quantum systems 155115611 . 
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SUPPLEMENTAL MATERIAL 



LINEAR QUANTUM SYSTEMS AND OPTIMAL 
UNRAVELLINGS 



In this section we recall the description of the evolution un- 
der continuous measurements for linear quantum system and 
we add some additional information regarding their optimal 
unravelling. 

We start by dealing with time-independent quadratic Hamil- 
tonians 



F = (l/2)R'^iJR , R=(£i,pi,, 



, Xn,Pn) 



As usual, the system is considered to be open and linearly 
coupled to a large uncorrelated Markovian bath. As explained 
in the paper, all the dynamics is defined by the 'drift matrix' 
A and by the diffusion matrix D, such that, in terms of the 
covariance matrix of the system, we have 



da- 



Act + a- A' 



D. 



(6) 



Given these matrices, and the Hamiltonian H, we can define 
the matrix 



M 



now 



i{H + nA), 



(7) 



and one can show, given the properties of A, D, and H, that 
this matrix is always hermitian and positive. Thus it can be 
diagonalized as M ~ VAV'', where V is unitary and A is a 
diagonal matrix and we can define the matrix C = A^^^V\ 
and the vector of operators c = CR. For instance, by consid- 
ering the two-mode case {n = 2) and D ~ {2N + l)l4, the 
matrix C and the vector c read 



/ x/iVi + 1 ^^/ThTl 



s/2 





V 



\ 



-i\/Ni u u 

VN2 + I WN2 + 1 

VN^ -iVN^ ] 



c= |al^/iV+l,a|^/iV,a2V^f + l,4^/^} ■ 



(8) 
(9) 



We now assume that the degrees of freedom of the environ- 
ment can be continuously monitored, resulting in instanta- 
neous Gaussian measurements on the system. These mea- 
surements, which are referred to as general-dyne [111, can be 
parametrized by an unravelling matrix C/, 



t/= - 



e + Rc[T] 
Im[T] 



Im[T] 
e - Rc[T] 



(10) 



where 8 and T are square matrices typically of dimension 
2n X 2n. To be a valid measurement, \J has to be posi- 
tive semi-definite, [/ > 0, and has to satisfy the following 
other constraints: T = T^ and = diag(?yi, . . . , 771,), where 



r]i G [0, 1]. The matrix U (or equivalently the pair of matri- 
ces O and T) corresponds physically to weak non-demolition 
measurements of the operators (c^S+c^T), which comprises 
all Gaussian measurements, including homodyne and hetero- 
dyne detection. In particular, the values i]i characterizing the 
matrix 9, correspond to the quantum efficiencies of the mea- 
surements on each environment mode. In our study we have 
always considered a common fixed efficiency rji = -q for each 
mode of the bath. 

Under such measurements, we obtain for the conditional state 
a diffusive equation with a stochastic component for the first 
moments (R)c and a deterministic equation for the CM <Jc, in 
formulae: 

d(R)c = ^(R)cdi+(crcC"^ + r"r)dw (11) 

^ ^Aa, + CT,A^ + D- {(T,C^ + T^)iCa, + T) 

(12) 

where dw is a vector of real Wiener increments satisfy- 
ing dwdwT ^ l2„df ii, r = (2C/)i/2S'C'fi, c-T = 
(Re[CT],Im[CT])and 



5* 



1„ 

-In 



Then a linear time-dependent term is added to the Hamiltonian 

Hf = -Il'nByit) , (13) 

where y(t) is the current obtained from the continuous mea- 
surement and the matrix B defines the kind of Markovian 
feedback action exerted on the system. In this case, the evolu- 
tion equation for the unconditional state covariance matrix cr 
is still of the form {d(T)/dt = A'cr + crA''^ + D' , where 

A' = A + BC , (14) 

D' ^D-C^B^ -BC + 2BB^, (15) 

andC ^2{Uy/^C . 

A CM (Tc is a physical stabilising solution of the conditional 
dynamics if it satisfies the two following conditions (see Eq. 
(fT2] i and notice that the second term on the right-hand side is 
always positive) 



An, + a^J^ + D>0 



(physicality condition) , (16) 
(stabilising condition) . (17) 



As derived by Wiseman and Doherty yj], given a stabilising 
CM (Tc, an optimal unravelling Uopt such that cTc can be ob- 
tained at steady-state, always exists. In particular a (not nec- 
essarily unique) optimal unravelling J7opt can be obtained by 
solving the equation 



2E'^UE = D + Ac 



.A\ 



(18) 



where E = Cctc + SCVl. The Hamiltonian term in Eq. (fTsT l 
is then chosen so as to cancel out the first moments and make 



the average unconditional state coincide with the conditional 
state. It can be shown that the matrix Sopt achieving this, for 
a given steady-state CM a^., reads 



-Sopt = ^(TcC — r . 



(19) 



BOUNDS ON MAXIMUM ACHIEVABLE SQUEEZING AND 
ENTANGLEMENT 



In the following we derive the bound on the Gaussian entan- 
glement and squeezing achievable by means of any feedback 
strategy where the control action consists in linear driving. We 
present our main findings as three lemmas leading to two final 
propositions. 

Lemma 1 (Bound on smallest symplectic eigenvalue) The 

smallest partially transposed symplectic eigenvalue V- of a 
generic CM a is bounded from below as follows 



v^ 



> AJAl , 



(20) 



A| and X2 being the two smallest eigenvalues of a. 



Proof Notice that this proof can be found in ||4|] . We will 
reproduce it here to make our work self-contained. 

Henceforth, \v) will stand for a unit vector in the phase 
space r and {v\ will be its dual under the Euclidean scalar 
product. Also, given a bipartition of the modes into the 'first' 
I and the 'last' 771 modes, let us define the matrix T, represent- 
ing partial transposition in phase-space, as T = 1® © ct®™, 
Gz being the z Pauli matrix. Hence, the partially transposed 
symplectic form is defined as f2 = T^T 

The squared symplectic eigenvalue iPl is the smallest 
eigenvalue of the matrix cr^/^fi^crOcri/^: 

i?_ = min(^;|c^l/2f^T^^^l/2|^^,_ 

\v) 

For each \v), one can define the unit vector \w) — 
£la'^l'^\v) I \J {v\(t\v), such that {v\a'^l'^\w) = (due to the 
antisymmetry of Jl) and 

iPl = niin{v\(T\v) {w\(t\w) > min (w|(T|w)(w|(T|t(;) ~ AjAa . 



Proof Note that the uncertain^ relation (fTST i is equivalent to 
the two following conditions j^]^: 



(T^/^n^ancr^^^ > 1 , and cr > 



(22) 



Foranyl?;) G Tone can define |z) = fl(T^''^\v) / y^ {v\cr\v), so 
that the Robertson Schrodinger Inequality (l22T i can be recast 
as 



{v\(T\v){z\cr\z)>l V|i))€r. 



(23) 



We will now denote by \vj) the eigenvectors corresponding to 
the increasingly ordered eigenvalues of cr: (t\vj) = AJ|wj). 
Let us consider a vector \v) belonging to the subspace, which 
we shall denote Tk, spanned by the k smallest eigenvectors 
of (T {|vj)}, for J < k. Clearly one has (u|cr|ii) < Xj,. The 
inequality ( |23] | then leads to 

xI{z\(t\z) > {v\(T\v){z\a\z) > 1 y\v)erk, 

which must be satisfied by all the vectors \z) belonging to the 
fc-dimensional linear subspace flVk (defined as the subspace 
spanned by the k orthogonal vectors fllvk))'- 

XI{z\(t\z) > 1 V|z) e nVk . 

By Poincare Inequality |6], a vector \z) must exist in ilTk for 
which {z\(t\z) < XJ., such that AJ!a|; > 1. q 

As an immediate corollary of Lemma 2, one obtains 



t\t 



AlA; 



> 



1 



AJA^ 



(24) 



Lemma 3 (Bound on eigenvalues of steady state CMs) 

Let (Tc be a conditional CM at steady state obtained under 
continuous Gaussian measurements, diffusion matrix D and 
a Hamiltonian matrix H. The product of the two largest 
eigenvalues Aj and Aj ofcTc is bounded as follows: 



x\xi < (^i±^ 



4 ajofg 



(25) 



vt 



l^> 



\v),\w) 



where {a A are the (strictly positive) eigenvalues of 
{—A — A^) in increasing order, while {S'^} are the (strictly 
positive) eigenvalues of D in decreasing order 



The last equality is easily verified once {v\(t^''^\w) = and 
(T > are enforced, and completes the proof. q 

Next, the uncertainty principle entails: 

Lemma 2 (Uncertainty relation for CMs' eigenvalues) Let 

{ XA and { XA be, respectively, the 2n increasingly-ordered 
and decreasingly-ordered eigenvalues of an n-mode CM a. 
Then one has: 



AJAJ > 1 for 1 < i < n. 



(21) 



Proof Given the condition ( fTTI ). and given the eigenvectors 
of (Tc, |A{) and jAj) corresponding to X\ and Aj, we have 

X\{X\\-[A^A^)\X{)<{X{\D\X{) (26) 

A^(A^|-(A + ^^)|A^)<(A^|i?|A^) (27) 



By defining A 
we have 



(j4+ A^) and multiplying the inequalities, 

AtA^(At|I|At)(A^|I|A^) < (At|i?|At)(A^|Z?|A^) (28) 



then 



\i\i / max(„i^2>=o(wi|-D|z;i)(i'2|-D|w2) 

■^1^2 — "^ ^^ 

min(„i|„2)=o(wi|A)|ui)(u2|A)|iJ2) 
{6\ + 4)2 



< 



4 Qfjaj 



where we use 



mm 

(■!;l|i;2)=0 



|wi|A|wi)(w2|^|w2> > a\a\ 



max {v^\D\v,){v2\D\v^) < ( ^^^ 



Further, and more generally, one has: 



(29) 
(30) 

(31) 

(32) 

D 



Proposition 1 (Maximal unconditional squeezing) Let cr 

be the CM of a steady-state achievable by continuous weak 
general-dyne measurements and linear driving in a system of 
bosonic modes subject to a drift matrix A and Gaussian white 
noise with a diffusion matrix D. The squeezing A| is bounded 
by 



Al> 



a 



,t 



4 



(33) 



Proof From Eq. (l2Tl l we obtain the relation A{ > 1/A|, 
where \\ (A|) is the largest (smallest) eigenvalue of a CM cr. 
By considering a conditional CM at steady state and following 
the same line of reasoning used in Lemma [3] we obtain the 
following inequality \\ < Si/a[, which yields the inequality: 



a1 > ^ > ^. 

At Sf 



(34) 



As explained before, the unconditional state g that we ob- 
tain from our dynamics is a statistical mixture (with Gaus- 
sian profile) of different conditional states gr having the same 
CM (Tc and different first moments r = (R)c, in formulae 
g = J dr p{r)gr. As a consequence, the unconditional CM 
reads cr ~ ctc + t where r > is the classical covariance 
matrix of the first moments' distribution p{r). Thus the lowest 
eigenvalue of a is lower bounded by the eigenvalue of (Tc and 
the bound above is valid for the unconditional state. It is worth 
to remember that, given an optimal CM cTc which is a physical 
stabilising solution of the conditional dynamics, the bound is 
tight, since we can always find a Markovian feedback strategy 
such that g — gr=o, that is such that the unconditional state 
has CM (Tc and zero first moments. q 

Proposition 2 (Maximal unconditional entanglement) Let 

g be the CM of a steady-state achievable by continuous weak 
general-dyne measurements and linear driving in a system of 
bosonic modes subject to a drift matrix A and Gaussian white 
noise with a diffusion matrix D. The logarithmic negativity 



Ej'/{g) A3/ of any 1 versus (n 
bipartition of g is bounded by 



1) modes or bisymmetric 



Ej^{g) < max 



0, l0g2 



Si + 4 



'^l'^2 / 



(35) 



Proof The chain of Inequalities ( l20b . ( l24b and ( |25] | leads to 

> :""",^ . (36) 



i^l 



(4 + 4)2 ' 



which, in turn, constrains the maximal logarithmic negativ- 
ity achievable for states g^. conditioned by Gaussian mea- 
surements having a CM ac- In fact, by using the formula 
£V = max[0, — log(t>„)] , we obtain. 



Ej^f{gr) < max 



0,log2 



"i"!. 



(37) 



On the other hand the unconditional (Gaussian) state reads 
g ^ J dr p(r)pr; this implies that g can be obtained from 
the Gaussian state pr=o (having CM (Tc and vanishing first 
moments) by local operations and classical communication 
alone, because first moments can be arbitrarily adjusted by 
local unitary operations. Since the log-negativity is an entan- 
glement monotone [8], we have E_\f{g) < E_\f{gr), that is the 
bound above is valid also for the unconditional state and can 
be achieved by means of optimal Markovian feedback. q 



Necessary conditions for the tightness of the bounds 

By considering how our bounds were derived, and working 
backward, we can determine sharp conditions on the matri- 
ces A and D, which fully characterise our dynamics, for the 
bounds to be achievable. In order to express such conditions, 
let us define the eigenvectors \a-) and |(5|) associated, respec- 
tively, to the j—th smallest eigenvalue of A = —A — A^ and 
j-th largest eigenvalue of D. This leads to the following two 
additional propositions: 

Proposition 3 (Conditions for maximal squeezing) A con- 
tinuously measured and linearly driven Gaussian system is 
able to saturate the bound ( 1331 ) if and only if 



J\ - 



I4> 



(38) 



Proof The inequality Aj < (5|/a| (the analogous of ( l30l l for 
the squeezing case) is only saturated if the eigenvector asso- 
ciated to the largest eigenvalue of cr coincides with \a[) and 
1 4)- hence our condition ( [38] l. in that it is always possible to 
construct a physical a with largest eigenvalue along a partic- 



ular direction. 



D 
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Proposition 4 (Conditions for maximal entanglement) A 

continuously measured and linearly driven Gaussian system 
is able to saturate the bound f |37| ) if and only if the following 
relationships are satisfied 



|4>T 



a. 



a. 



V2 
\5\) ± |<5: 



{al\T\al\ 



V2 
= 



(39) 

(40) 

(41) 
(42) 



(where =p and ± mean f/zaf if Eq. ( 1591 ) /zoi a minus sign then 
Eq. ( 1401 ) /jfli fl pZwi, fl«rf v/ceverifl, fl«t/ that either choice is a 
valid condition). 

Proof Eqs. (|39] l and ( l40b are necessary for the saturation 
of Inequality ( l30b , along with the choices |A{) = |q;|) and 



k = 1, k = 2 and, by induction, for any fc, reveals that the 



condition A|;A|^ 



1 is saturated if and only if |A^) = 17 1 A 



i\ 



it 



(where the eigenvectors associated to Xj, and A|, have been 
denoted with \\j,) and |AJ.)). Now, in order to saturate the 
bound, this additional condition can only be imposed if the 
two eigenvectors |A|;), already determined by (|39t and (|40] |. 



are orthogonal to i7|AJ;) 
vectors of cr), that is 



(so that the latter can also be eigen- 



(A^I^IAt) = 



1^1 



= 



(43) 



Further, inspection of Lemma 1 shows that, for Inequality (l20l i 
to be saturated, it must be jAj) ~ f^|A|) which, by the con- 
ditions I At) = |a|), |A^) == |a^) and |A^) = n\Xi) imposed 
at previous steps, becomes lal) — fl^flfllal), which proves 
condition (l4Tl i. Finally, by inserting Eq. (l4Tl i into Eq. ( |43] |, 
and noting that il^fiilil = — T, one can recast condition (l43l l 
in terms of |a|) alone as condition (|42] |. q 



APPLICATIONS. 
OPTIMAL ACHIEVABLE STATES AND UNRAVELLINGS 

In this section we give some details regarding the optimal 
states and the optimal unravellings in the different cases we 
discuss in the paper. 

Free systems. - Let us start by considering the maximum 
single-mode squeezing achievable when no Hamiltonian is 
present. For a given number of thermal excitations N, one can 
achieve a squeezed vacuum state \ip) = S'(ropte*'^)|0), where 



5(e)=cxp{ie(an' 
operator and 



i^*a^} is the single-mode squeezing 



ropt = iln(^ 



2N 



(44) 



Since there is no privileged direction in the phase-space, one 
can choose the phase of the squeezing cf) by choosing one 
of the optimal measurements on the environment. In fact 
the optimal unravelling and feedback strategy correspond to 
perform a homodyne measurement on a quadrature x^ = 
cos (f> X + sin (f) p, and then to drive the orthogonal quadra- 
ture with the obtained current. For instance, if one measures 
X, the optimal unravelling matrix reads 



C/o 



ipt 



/O 



Vo 






N{l + N) 
1+N+N'-' 









N{l+N) 

1 + N+N 



\ 



u 



(45) 



and the steady-state will be squeezed in the x quadrature. The 
two matrices that parametrize the measurement, are 



e. 



jofj). Then, inspection of the proof of Lemma 2 for where 



opt 



'7opt 



-T 



opt 



jyoptls 



N{l + N) 
1 + N + N^' 



(46) 



(47) 



As explained in the paper, this corresponds to an inefficient 
homodyne measurement, with efficiency r/opt- 
As regards the entanglement, we focus here on the two-mode 
case where the two baths have the same number of thermal 
excitations N. We find that a two-mode squeezed state |i/;) ~ 
'S'2(?'opte"^)|0), is a stabilising solution of the dynamics and 
achieves maximum value of entanglement given by the bound 
we have derived. We denote with S-ziO ^ exp{^a^^6^ — ^*a5} 
the two-mode squeezing operator, while the squeezing param- 
eter Topt is still equal to the one given in (l44l) . As in the pre- 
vious case, there is no privileged direction in the phase-space 
and thus one can consider different optimal unravellings, giv- 
ing different squeezing phases 0. One of them corresponds to 
the measure of the quadratures x\ — X2 and p\ +p2, which in 
terms of the unravelling matrix reads 



C^opt - 2 



/ 1 



-1 





V 





1 



-1 








-1 



1 









o\ 
-10 


10 10 
10 1 



10 10 

1 ly 



(48) 



By looking at the corresponding matrices 



Topt = 



e = I4 



(49) 
(50) 



one observes that in this case the optimal unravelling requires 
a perfect homodyne detection on the non-local quadratures. In 
the results we present in the paper, we also consider the role of 
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the efficiency on the achievable entanglement; this can be in- 
vestigated by just replacing the optimal unravelling, with the 
matrix U = rjUopt, where rj G [0, 1] denotes the homodyne 
efficiency which is considered to be equal for both modes in- 
volved. 

Parametric Hamiltonians. - We consider here the case of 
degenerate parametric down conversion, described by the 
quadratic Hamiltonian 

between two modes at the same frequency, such that the av- 
erage number of thermal excitations in the two-modes are the 
same and set equal to N . In this case the coupling defines 
a privileged direction in the phase space, since it squeezes 
the quadratures xi — £'2 and pi + p2- Indeed, the two- 
mode squeezed state which is stabilising solution of the dy- 
namics and achieves the bound on the entanglement, reads 
<S'2(n)pt)|0), where 



'opt 



1 fl + 2N 
— In 



(51) 



2 V 1 - 2x . 
The state is squeezed, as expected, in the two non-local 



quadratures indicated above. Also the optimal unravelling that 
we obtain with Eq. ( fTSl ) reflects this fact. It corresponds to the 
matrix given in Eq. ( |48] |. and thus to a perfect homodyne mea- 
surement of the quadratures xi — X2 and pi +p2- In this case 
too, results for inefficient measurements, that is by consider- 
ing the unravelling U = rjUopt, are reported in the Letter 
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